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A TRANSITION MATRIX APPROACH 
TO THE DAVENPORT GYRO CALIBRATION SCHEME* 

G. A. NatansorV 

The in-flight gyro calibration scheme commonly used by NASA Goddard Space 
Flight Center (GSFC) attitude ground support teams closely follows an original 
version of the Davenport algorithm developed in the late seventies. Its basic 
idea is to minimize the least- squares differences between attitudes gyro- 
propagated over the course of a maneuver and those determined using post- 
manuever sensor measurements. The paper represents the scheme in a recursive 
form by combining necessary partials into a rectangular matrix, which is 
propagated in exactly the same way as a Kalman filters square transition matrix. 

The nontrivial structure of the propagation matrix arises from the fact that 
attitude errors are not included in the state vector, and therefore their derivatives 
with respect to estimated gyro parameters do not appear in the transition matrix 
defined in the conventional way. 

In cases when the required accuracy can be achieved by a single iteration, 
representation of the Davenport gyro calibration scheme in a recursive form 
allows one to discard each gyro measurement immediately after it was used to 
propagate the attitude and state transition matrix. Another advantage of the new 
approach is that it utilizes the same expression for the error sensitivity matrix as 
that used by the Kalman filter. As a result the suggested modification of the 
Davenport algorithm made it possible to reuse software modules implemented in 
the Kalman filter estimator, where both attitude errors and gyro calibration 
parameters are included in the state vector. 

The new approach has been implemented in the ground calibration utilities used 
to support the Tropical Rainfall Measuring Mission (TRMM). The paper 
analyzes some preliminary results of gyro calibration performed by the TRMM 
ground attitude support team. It is demonstrated that an effect of the second 
iteration on estimated values of calibration parameters is negligibly small, and 
therefore there is no need to store processed gyro data. This opens a promising 
opportunity for onboard implementation of the suggested recursive procedure 
by combining it with the Kalman filter used to obtain necessary attitude 
solutions at the beginning and end of each maneuver. 


* This work was supported by the National Aeronautics and Space Administration (NASA) / Goddard Space Flight 
Center (GSFC), Greenbelt, Maryland, Contract GS-35F-4381G, Task Order No. S-03365-Y. 
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I. INTRODUCTION 


Propagation of a state vector front one measurement time to another is usually done 1 by 
introducing a transition matrix formed by partial derivatives of the current state with respect to 
the state at an epoch time. A well-known technique 2 - 3 has been developed to propagate the 
transition matrix between sequential measurements using gyro data. The paper extends this 
propagation technique to the Davenport gyro calibration scheme . 4-7 The main obstacle to such 
an extension comes from the fact that the cited gyro calibration scheme treats an a priori given 
change in the spacecraft attitude within a specified time interval as a pseudo-measurement, and 
therefore attitude errors are not included in the state vector anymore; as a result, their derivatives 
with respect to estimated gyro parameters (such as misalignments, biases, and scale factors) do 
not appear in the transition matrix defined in the conventional way. To overcome this 
complication, the new approach combines necessary partials into a rectangular matrix, which 
can be propagated in exactly the same way as the conventional ( square by definition) transition 
matrix . 2 - 3 

Assuming that the first iteration eliminates bulk errors, representation of the least-squares gyro 
calibration scheme in a recursive form allows one to discard each gyro measurement 
immediately after it is used to propagate the attitude and state transition matrix. Due to a 
significant decrease in required storage size, this feature of the new approach seems especially 
promising for onboard applications. 

The next Section presents a simplified derivation of the original version of the Davenport 
algorithm . 4 - 5 Its final result is an explicit expression of the vector attitude residual in terms of the 
error sensitivity matrices iy k utilized by Kalman filter estimator . 3 It is shown that the derived 
expression turns into the conventional one 7 if only linear terms are kept in the expansion of each 
matrix y k as a Taylor series in the duration At k of the k ^ 1 propagation interval. 

Section III introduces a rectangular matrix which gyro-governed evolution is performed via the 
same recurrence relations as those used for gyro propagation of the conventional state transition 
matrix. The derivation is accomplished in Section IV, which outlines main steps of the 
suggested recursive procedure. 

The new algorithm has been implemented* and successfully used to calibrate gyros 9 for the 
Tropical Rainfall Measuring Mission (TRMM). One of the advantages of the suggested 
modification of the Davenport algorithm is that it allowed a reuse of software modules 
implemented in the Kalman filter estimator , 3 where both attitude errors and gyro calibration 
parameters are included in the solve-for state vector. Section V discusses some preliminary 
results of the TRMM gyro calibration. It is shown that the required accuracy of gyro calibration 
can be indeed achieved by a single iteration, and therefore each gyro measurement can be indeed 
discarded immediately after it was used. The paper also studies a possibility to reduce a volume 
of processed gyro data without jeopardizing the accuracy of calibration. 
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II. MATHEMATICAL GROUNDS OF DAVENPORT ALGORITHM 


The Davenport method 4 is a two-step procedure. The first step is to determine attitudes at the 
ends of specially selected calibration intervals. For successful calibration, the selected time 
intervals usually cover a series of maneuvers associated with significant changes in body rates. 
It is essential that, regardless of maneuver specifics, each calibration interval must both start and 
end in a constant-rate mode. To determine the spacecraft attitude at the ends of each interval, 
sensor measurements are then collected only during time periods within constant- rate modes, 
when unknown errors in gyro misalignments and scale factors are compensated by additional 
gyro biases estimated simultaneously with the spacecraft attitude. As a result, one can assume 
that gyro propagation from one sensor measurement time to another is done accurately enough, 
despite the fact that gyros have not been properly calibrated yet. 

The second step is quaternion propagation starting from the predetermined attitude quaternion at 
the beginning of each calibration interval and stopping at its end. The resultant propagated 
quaternion is then compared with the second of two attitude quaternions predetermined for this 
calibration interval. The comparison is done by multiplying one of the two quaternions at the 
ending time by the inverse of other. The vector part of the product is then treated as a vector 
residual, with the total number of these attitude quaternion vector residuals (AQVRs) always 
equal to the number of the calibration intervals. 

Davenport’s principal result is an approximate expression for the AQVRs in terms of vector 
deviations of the observed angular velocity vectors from the true rates. The outline of the 
Davenport method presented here mainly follows Keat’s 5 interpretation of Davenport’s original 
work. 4 To simplify the notation, the discussion will be limited only to a single maneuver so that 
the index labeling different maneuvers can be omitted. An extension of the final expression for 
an AQVR to a series of sequential maneuvers is performed in a trivial way by attaching an 
additional index to both residual and all angular velocity vectors. 

Let be an observed angular velocity vector obtained by adjusting measured rates with some 
estimated parameters, where subscript k refers to the k-th available gyro measurement within the 
maneuver in question. The observed vector differs from a true vector CO ^ by a rate error 5co ^ , 
that is. 


co k = cojf j - 665 k . 0) 

Both vectors <3 and (5 k are assumed to remain constant during a time interval At k so that the 
quaternion propagated over n intervals (starting from the known quaternion ) can be 
represented as 

Vp = q M n (2) 

k=l 
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where 10 q((5 At) = |cd sin(<J>/2), cos(<|>/2)J, with to = co / [63 1 and <j> = [65 1 At . The AQVR 
Z is defined via the relation:* 

[Z.-fiW ( 3 ) 

where q fin is the given attitude quaternion in the end of the maneuver. 

It is assumed that the attitude quaternion q fin can be obtained from q jnit by propagating the 
latter with the true constant angular velocity vectors CD k over the time intervals At k , so that 

flfin ~ *Iinit I! Q(® k At k ) . (4) 

k=l 


To express the AQVR Z in terms of errors in gyro parameters, one first needs to linearize the 

quaternion product q f jj, q prop in 5(0 k . At this point one has to deal with unnormalized 

quaternions, which form the so-called ‘associative algebra’. Note that both attitude and 
propagation quaternions discussed above are normalized quaternions, which cannot be either 
summed up or multiplied by a scalar, in contrast with unnormalized quaternions. On the other 

hand, the inverse operation q" 1 is well defined just for normalized quaternions. Only the 
multiplication law given by Eq. (D-8) in Ref. 10 is common for both normalized and 
unnormalized quaternions. It is essential that, by analogy with orthogonal matrices, the 
multiplication law is associative, i. e., q(q'q") = (qq')q" for any three unnormalized 
quaternions, q,q\ and q". Another important features of the multiplication law are that 
q(q' + q") =qq' + qq" and that q(kq') =(kq )q' for a scalar multiplication. After the 
mentioned features of the multiplication law are established, unnormalized quaternions can be 
formally treated in the same way as square matrices, with the norm of quaternion given by Eq. 
(D-9) in Ref. 10 used instead of the matrix determinant. In particular, all Taylor expansions look 
very similarly, except that each product should be computed based on the quaternion 
multiplication law. 

Substituting Eqs. (2) and (4) into Eq. (3) and keeping only terms linear in 5© k , one can 
represent the latter expression as 

8Q*(l-n)[0,ll+£ q,!.„8q k q^„ (5) 

1 J k=l 

where 

k-f-Ak _ 

Qk^k+Ak s n <j(G> k- J At k .) for Ak =1,..., n-k > 0, (6a) 

k-k+l 


* Note that our definition of the AQVR differs by the factor (-1) from that used by Keat. 5 
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and 


(6b) 


q„->n =[0> !]. 

ik.VH^I 2 ] = 5qk - q(-® k At k)q(® k dJAt k) • ( ? ) 

Note that the sum in the right-hand side of Eq. (5) is formed by the products of normalized 
quaternions, and hence, to simplify each product, one can take advantage of the existent 
isomorphism between normalized quaternions and 3x3 orthogonal matrices. Making use of 
Eq. (12-7b) in Ref. 10, one can easily verify that 

RA(g)R T = A(Rg), (8) 

where R is an arbitrary 3x 3 orthogonal matrix, whereas the rotation matrix A(g) associated 
with the Gibbs vector g is given by Eq. (12-7b) in Ref. 10, with g= etan(<D/2). Representing 
Eq. (8) in the quaternion form and substituting the resultant expression in the sum in the right- 

hand side of Eq. (5), one can finally represent the AQVR Z as 

Z=z [R^„i k -0(|5(5»| 2 )], (9) 

where R^k 1 ' s the rotation matrix associated with the propagation quaternion q k _, k . 

As discussed in detail in Ref. 1 1, an explicitly expression of the vector z k in terms of the rate 
error 5 CO k has the form: 

z k = 5c5 k + 0(|5(5 k | ), (1®) 

where V[/ k is the error sensitivity matrix used by the Kalman filter estimator, 3 that is, 

vj/ k =I 3 At k -^[6 a k dj x]At 2 k v 2 (^/2) + 2[o5 a k dj x] 2 ri((j^j)At^ (11) 

with v((p) s sin <p/cp , = 05 (1 — sincp)/ cp 2 , and s|a>|^|At k . Note that a slightly 

different representation for error sensitivity matrix (11), compared with Ref. 3, makes it possible 

to compute this matrix using expansions explicitly stable at the limit Icoj^l — * ► 0- 

Finally, AQVRs (9) are represented as linear combinations of errors Ax, (i=l,...,p) in gyro 
parameters: 

Z» < 12 > 

2 i=t k=i 
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(13) 


by substituting the rate errors 

5d5 k = £ fjfwfjAXi . 

into Eq. ( 1 0). Computation of necessary partials is then performed in a trivial way. 

Note that Keat’s formula 5 for the AQVR (utilized in the conventional version 7 of the Davenport 
algorithm 4 ) is obtained from Eq. (13) by keeping only the first term in the right-hand side of Eq. 
(11), which seems to be a sufficiently accurate approximation in most cases (see comments made 
in the end of Section V). Another minor modification comes from a slightly different choice of 
the state vector Ax , which is formed by three bias errors Abj (i=l,2,3), three scale factor errors 
Akj (i=l,2,3), and two misalignment angle errors e k j (k=l,2) for each of three gyros (i=l,2,3). 

Such a choice of gyro calibration parameters* makes it possible to calibrate each gyro 
separately, which is convenient in case of a spacecraft having only one gyro, such as Solar and 
Heliospheric Observatory (SOHO). 12 


III. PROPAGATION OF ATTITUDE MATRIX VECTOR RESIDUALS 

The main purpose of this Section is to show that an attitude residual can be represented in the 
general form: 


0 = HO staIe Ax , (14) 

where 0 , Ax , H , and O sme are usually referred to as a measurement residual, a state error 
vector, a sensitivity matrix, and a state transition matrix, respectively. The crucial point is that 

the transition matrix Q c „„ can be computed as the last term O state = O n in a sequence of 
recurrence relations: 

= ^k-l-^k 3>k-i * ( 1 5 ) 

where the (p +3) x (p+3) matrix O k _,_, k is an incremental transition matrix conventionally used 
in Kalman filter applications to propagate the attitude state vector (see, for example. Eq. (F8-26) 
in Ref. 3). A certain complication, however, comes from the fact that the calibration scheme in 
question estimates only gyro calibration parameters, and therefore attitude errors are not 

included into the state vector Ax . As a result O k turn out to be rectangular matrices having 
p+3 rows but only p columns. 

As mentioned above, the initial and final attitudes in the Davenport method are determined using 
sensor measurements in constant-rate modes, when solved-for biases compensate for errors in 

gyro misalignments and scale factors. The attitude matrix vector residual (AMVR), 0 , is then 
defined via the relation: 


378 



exp[0 x] = A obs A prop , 


(16) 


where A obs and A prop are the attitude matrices associated with the attitude quaternions q fin 
and q prop in the previous Section.* One thus finds 

0 sin(|0|/2)/|0|=Z . (17) 

Since we are interested only in terms linear in gyro errors, one can simply put 

0 = 2Z (18) 


and make use of Eq. (12) to represent the AMVR 0 as the last term0 ~9 n in the sequence 


i=l v ' 


(19) 


with 0 O =0 . 

Eq. (19) immediately leads to the conventional Kalman filter expression 3 for propagation of the 
combined attitude error / gyro calibration parameter state vector: 


0i 


Ax 


= <D, 


k-l-»k 


©k-1 

Ax 


where 


R 


0 ), 


k-l-»k 


¥k F k 


k-l->k 


0 


px3 


with 




By initializing sequence (15) via the relation 


( 20 ) 


( 21 ) 


( 22 ) 


* The author is thankful to J. Sedlak for pointing to a misprint in the definition of the AMVR0 in Ref. 1 1 
leading to a sign error in Eqs. (III-5), (III-6), (III-8), (HI-14), and (III-17) there. 
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(23) 




one finds that 



V1F1 

! P 


(24) 


and hence, making use of Eq. (20) at k=l, 

0! = HO ,Ax, (25) 

with 

H-[l 3 «3x P ]- (26) 

By applying mathematical induction to Eq. (20) and making use of the fact that the first three 
rows of the matrix HO k coincide with the first three rows of the transition matrix O k for any 
k, one can easily verify that 

§ k =HO k Ax. (27) 


Substituting the state transition matrix O state for O n then immediately leads to Eq. (14), which 
constitutes the main result of this work. 


IV. REPRESENTATION OF THE DAVENPORT GYRO CALIBRATION SCHEME IN A 
RECURSIVE FORM 

The suggested recursive procedure has been implemented in the following way.* 

Estimation starts by setting elements of the so-called ‘measurement accumulation’ vector Au to 
zero. One also initializes elements of the covariance matrix P with some a priori values. The 
attitude matrix is then propagated from one gyro measurement to another: 

■^k — ^k-l-»k^k-l (28) 

starting from the given observed attitude A 0 associated with the quaternion q init at the 
beginning of the first maneuver. At the same step one also computes error sensitivity matrix 
(11) and rate-dependent partials in the right-hand side of Eq. (22) which are then substituted. 
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together with the rotation matrix R^j. , into Eq. (21) for the matrix used to propagate 

the state transition matrix via recurrence sequence (15). At the end of the maneuver the AMVR 

0 is computed by linearizing Eqs. (16): 

0 ~ — [5A23 - 5Aj2 , SAji — 6Aj3 , 6A,2 5A 2 i] 

(29) 

where 5Ay are elements of the orthogonal matrix 


= ^obs ^prop 

(30) 

After the AMVR 0 is computed, one updates the measurement accumulation vector and the 
inverse covariance matrix W = P* 1 according to the standard equations: 

Au <— Au +O s ^ te H T C _1 0 
and 

(31) 

W W + O s { aIe H T C" 1 HO state , 

(32) 


where C is a 3x3 measurement covariance matrix. The state transition matrix is then reset to 
O 0 and the processing continues starting from the beginning of the next maneuver. 

After the last maneuver is processed, one obtains the covariance matrix P by inverting the 
resultant W matrix and computing the state error vector Ax from the measurement 
accumulation vector Au : 


Ax=PAu. (33) 

The magnitude of the state error vector Ax is then compared with the given tolerance to proceed 
with iterations if necessary. 


V. TRMM IN-FLIGHT GYRO CALIBRATION 

The TRMM is an Earth-pointing three-axis stabilized spacecraft. Its body z-axis is nominally 
pointed along the geodetic nadir. 13 It can be in ‘+x forward’ or ‘-x forward’ nominal mode, with 
its body x-axis being approximately either parallel or anti-parallel to the spacecraft velocity. For 
power and thermal protection of science instruments from direct exposure to the Sun, yaw 
maneuvers from one nominal mode to another (similar to those depicted in Figs, la and lb) are 
periodically performed. Since the body y-axis is parallel (anti-parallel) to the orbit normal in the 
‘-x forward’ (‘+x forward’) mode, the only nonzero component of the spacecraft angular velocity 
vector is the pitch rate, nominally equal to +1 revolution per orbit (RPO) in the -x forward 
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mode and -1 RPO in ‘+x forward’ mode. This is illustrated by the flat portions of the dashed 
curves in Figures la and lb. It is essential that rates remain nearly constant in both nominal 
modes, and therefore the spacecraft attitude can be determined with a sufficient accuracy without 
a complete gyro calibration, provided that gyro biases are included in the state vector to be 
solved for. (Unknown errors in scale factors and gyro misalignments manifest themselves as 
some additional biases, which differ for different modes.) 



Time (sec) 


Figure 1 . TRMM body rates during +X to —X (upper) and —X to +X (lower) yaw maneuvers. 

On December 14, 1997 the TRMM was also placed in the *-y forward’ mode to calibrate 
scientific instruments, with the y body axis being anti-parallel to the spacecraft velocity vector. 
As seen from gyro rate profiles depicted in Fig. 2, the spacecraft was rotating for about one hour 
around its body x-axis with the rate of +1 RPO, before coming back to the ‘-x forward’ mode. 
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In addition to maneuvers between the three Earth-pointing modes mentioned above, the TRMM 
was commanded on January 7, 1998 to stay for one orbit in the inertial hold mode used to 
calibrate a science instrument by pointing it toward cold space. Fig. 3 presents the 
corresponding x and y body rates. (The z body rate is omitted since its deviations from zero 
would be practically invisible at the figure scale.) 



Table 1 lists time intervals selected for calibration of the TRMM inertial reference unit (IRU).’ 
The reference attitudes were determined at the beginning and at the end of each maneuver using 
Digital Sun Sensor (DSS) and Barnes Static Earth Sensor Assembly (SESA) measurements.' 3 

Table 1 - Intervals of Gyro Data Used for IRU Calibration 


Maneuvers 

Time Interval 

Spacecraft Activity 

1 

9:45:00- 11:25:00, Dec. 13, 1997 

+x forward 

2 

13:11:00- 13:29:00, Dec. 13, 1997 

+x to -x yaw maneuver 

3 

15:05:00- 16:40:00, Dec. 13, 1997 

-x forward 

4 

11:43:30- 12:58:00, Dec. 14, 1997 

-y forward mode 

5 

13:07:00 - 14:39:15, Jan. 7, 1998 

Inertial hold 

6 

20:47:30 - 21:06:00, Jan. 14, 1998 

-x to +x yaw maneuver 


Three residuals per maneuver were then obtained by propagating the spacecraft attitude with 
gyro rates from the beginning of each maneuver and comparing the result with the predetermined 
reference attitude at the end of the maneuver. Table 2 presents roll, pitch, and yaw attitude 
residuals obtained by gyro propagation with pre-launch and calibrated gyro biases, scale factors, 
and misalignment matrix. The corresponding values of calibration parameters for each gyro 
(i. e., a bias, a scale factor, and misalignment angles* relative to body axes) are listed in Table 3. 

" A deviation of the gyro axis from its nominal direction can be derived from two other misalignments; 
however, it is included in Table 2 just to simplify the notation, with zeros standing for some rather small 
numbers completely irrelevant to our discussion. 
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Table 2- Roll, Pitch, and Yaw Attitude Residuals (deg) Before and After Calibration 


Maneuvers 

At = 0.5 s 

Pre- launch 
At = 1 s 

At = 2 s 

At = 0.5 s 

1 st iteration 
At = 1 s 

At = 2 s 

i 


m mm 


■ 












H 


0.002 



2 

0.195 

0.212 

0.202 





0.002 

-0.001 

-0.003 



-0.001 


0.024 

0.031 

0.031 



0.002 

3 

0.018 

0.015 

-0.021 



mxm 


-0.433 

-0.412 

-0.441 



-0.003 


-0.011 

0.003 

0.034 



-0.042 

4 

0.235 

0.233 

0.265 


0.001 

-0.001 


-0.243 

-0.233 

-0.243 

1 

-0.001 

-0.001 


-0.108 

-0.116 

-0.131 


-0.006 

0.007 

5 



0.099 

mMm 


-0.001 


-0.194 


-0.248 



0.011 




0.075 



0.006 

6 



0.271 



0.006 



0.086 

0.078 



0.003 


0.047 

0.045 

0.074 



0.004 


Table 3 - Calibration parameters 



Pre-launch 


1st iteration { 



At = 0.5 s 

& 

ii 

(A 

At = 2 s | 

Biases 





(deg/sec) 


1 


1 


WmMwM 



1 

Scale factors 

1.00000 

1.00044 

1.00035 

1.00031 


1.00000 

1.00053 

1.00046 

1.00046 


1.00000 

1.00088 

1.00092 

1.00096 

Misalignment angles (deg) 

0.000 

0.000 



relative body axis x 

0.059 

0.044 


1 


0.005 

0.054 


| 

Misalignment angles (deg) 

-0.036 




relative body axis y 

0.000 


0.000 



0.030 




Misalignment angles (deg) 





relative body axis z 

I 

-0.124 




| 
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To investigate the possibility of reducing the amount of processed gyro data (about 60 000 points 
for all six maneuvers), propagation was performed using several time steps: At = 0.5 s, 1 s, 2s. 
Inspection of Table 2 and Table 4 shows that every other gyro measurement can be successfully 
skipped without any noticeable effect on the accuracy of estimation. Skipping three of every 
four gyro measurements still gives reasonably good results, though some degradation in the 
accuracy can be clearly seen. 

Table 4 - Mean and RMS Residuals (deg) 



Pre- launch 

At = 0.5 s At = 1 s At = 2 s 

1st iteration 

At = 0.5 s At = 1 s At = 2 s 

Mean 

Residuals 

0.1343 0.1404 0.1343 

-0.1343 -0.1426 -0.1603 

0.0103 0.0105 0.0226 

■ 

RMS 

Residuals 


I 


It is essential that the calibration be accomplished by a single iteration. This can be easily seen 
by comparing corrections to gyro biases and elements of the G-matrix due to the first and second 
iterations, presented in Table 5. Contributions to attitude residuals from the second iteration are 
so small that they would have no effect on the values t in Tables 2 and 4 (to the precision 
shown). 


Table 5- Bias and G-Matrix Corrections 



1st iteration 

x-gyro y-gyro z-gyro 

2nd iteration 

x-gyro y-gyro z-gyro 

Bias 

corrections 

(deg/sec) 

-1.74 xlO" 5 4.15 xlO* 5 -1.15 xlO* 5 

-3.89x10*® 4.18x10® 0.11 xlO* 8 

Corrections 
to G-matrix 
elements 

-0.44 xlO* 3 0.39 xlO* 3 -0.57 xlO" 3 
0.27 xlO" 3 -0.54 xlO* 3 1.59 xlO* 3 
-0.86 xlO' 3 -1.16 xlO" 3 -0.88 xlO* 3 

-2.69 xlO* 6 -1.90 xlO"* -0.72 xlO* 6 

-2.86 xlO* 6 -0.58 xlO"* 0.81 xlO* 6 

-0.85x10"* -4.08x10"* -1.43x10"* 


To study the significance of higher-order powers of At* in the error sensitivity, calibration was 
repeated using only the linear terms in Eq. (11), which is equivalent to the use of the Davenport 
method in its conventional implementation. 7 It was found that neglecting higher- order terms 
does not practically affect the calibration results, so that the linear approximation seems to be 
sufficiently accurate for calibration purposes. 
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VI. CONCLUSIONS 


A new approach has been implemented in the ground calibration utilities* and successfully used 
for the TRMM IRU calibration. It has been shown that the required accuracy can be achieved 
by a single iteration. As a result the new approach seems to be especially useful for onboard 
applications by allowing one to discard each gyro measurement immediately after it is used to 
update the state vector and covariance matrix. 

Recently the Rossi X-Ray Timing Explorer (RXTE) ground launch support team has reported 14 
some problems in Kalman filter estimation of gyro scale factors and gyro misalignments, and the 
new least-squares approach to gyro calibration makes it possible to extend advantages of the 
Davenport algorithm to onboard applications. To avoid memory-consuming batch attitude 
determination, one can use the Kalman filter to determine the spacecraft attitudes before and 
after each of the selected maneuvers. The Kalman filter state vector is composed only of attitude 
and gyro bias errors, so that gyro biases will change with each new gyro measurement, in 
contrast with those used in the Davenport method. For this reason one has to propagate in 
parallel two separate transition matrices, namely, propagation rates for the spacecraft attitude and 
the transition matrix used by the suggested recursive algorithm are obtained by adjusting raw 
measurements with a priori biases which remain the same for all the selected maneuvers. On the 
other hand, propagation rates for attitude errors and the transition matrix utilized by the Kalman 
filter are obtained by adjusting raw measurements with the solved-for biases (and the same scale 
factors and misalignments as in the former case). Feasibility of this approach is currently 
investigated using Submillimeter Wave Astronomy Satellite (SWAS) simulated data. 
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